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Abstract 

Several studies suggest that one possible cause of impaired wound healing 
is failed or insufficient lymphangiogenesis, that is the formation of new lym¬ 
phatic capillaries. Although many mathematical models have been developed 
to describe the formation of blood capillaries (angiogenesis) very few have 
been proposed for the regeneration of the lymphatic network. Moreover, lym¬ 
phangiogenesis is markedly distinct from angiogenesis, occurring at different 
times and in a different manner. Here a model of five ordinary differen¬ 
tial equations is presented to describe the formation of lymphatic capillaries 
following a skin wound. The variables represent different cell densities and 
growth factor concentrations, and where possible the parameters are esti¬ 
mated from experimental and clinical data. The system is then solved nu¬ 
merically and the results are compared with the available biological literature. 
Finally, a parameter sensitivity analysis of the model is taken as a starting 
point for suggesting new therapeutic approaches targeting the enhancement 
of lymphangiogenesis in diabetic wounds. The work provides a deeper under¬ 
standing of the phenomenon in question, clarifying the main factors involved. 
In particular, the balance between TGF-/3 and VEGF levels, rather than their 
absolute values, is identified as crucial to effective lymphangiogenesis. In ad¬ 
dition, the results indicate lowering the macrophage-mediated activation of 
TGF-/3 and increasing the basal lymphatic endothelial cell growth rate, in¬ 
ter alia^ as potential treatments. It is hoped the findings of this paper may 
be considered in the development of future experiments investigating novel 
lymphangiogenic therapies. 
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1. Introduction 

1.1. Motivation 

Much effort has been spent in order to better understand and poten¬ 
tially treat the impairment of wound healing in diabetic patients. In this 
regard, one phenomenon that has recently gained attention from biologists 
is lymphangiogenesis] that is, the formation or reformation of lymphatic vas¬ 
culature m la El sj. Insufficient lyphangiogenesis, as observed in diabetic 
subjects, appears to correlate with failed or delayed wound healing. 

Impaired wound healing is a major health problem worldwide and in 
recent decades has attracted the attention of both biologists and mathemati¬ 
cians. In many cases unresolved wound healing correlates with prolonged 
infection, which negatively affects the patient’s quality of life, causing pain 
and impairing their physical abilities. Particularly serious infection may even 
require the amputation of an affected limb [5]. Furthermore, impaired wound 
healing also constitutes a major problem for health care systems, accounting 
for approximately 3% of all health service expenses in the UK and 20 billion 
dollars annually in the USA laiaiHii. Several systemic factors contribute to 
the delay or complete failure of the wound healing process [iDiiiDiiaiiaE!. 
In particular, diabetic patients exhibit a slower and sometimes insufficient 
response to infection after injury. Such a delay often results in a chronic 
wound; that is, the wound fails to progress through the normal stages of 
healing and usually remains at the inflammation stage [laiin]. 

Interest in lymphangiogenesis in reference to wound healing is very recent: 
for example in the Singer & Clark 1999 review m the process is not men¬ 
tioned. Nevertheless, lymphatic vessels have recently become regarded as a 
crucial factor in wound healing [iHiiin!. They mediate the immune response 
and maintain the right pressure in the tissues 120 !, thus playing a very impor¬ 
tant role in inflammation and contributing to the healing of a wound [211 [22]. 
Moreover, failed restoration of a lymphatic network (observed, for example, 
in diabetic patients) is now thought to be a major cause of impairment to 
wound healing uni ESI ED- 

Mathematical modelling has proven a useful tool in understanding the 
mechanisms behind numerous biological processes. It is therefore of interest, 
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and potentially great utility, to build a model describing lymphangiogenesis 
in wound healing, considering both the normal and pathological (diabetic) 
cases. 

1.2. Biology 

Wound healing is a very complex process involving a number of entwined 
events, which partially overlap in time and influence one another. For sim¬ 
plicity and educational purposes, it is often divided into four different phases: 
hemostasis, inflammation, proliferation and remodeling. Here the key events 
in each of the phases are summarised; for further details, see for instance 
0 ESI EBl [m El!). A few minutes after injury, the contact between blood 
and the extracellular matrix (ECM) causes a biochemical reaction that leads 
to the formation of a blood clot. This “crust” has the double function of 
stopping the bleeding (hemostasis) and providing a “scaffold” for other cells 
involved in the process to be described below. Concurrently, chemical reg¬ 
ulators (such as Transforming Growth Factor /3, or TGF-^) are released, 
which attract cells such as neutrophils and monocytes to the wound site. 
These cells clean the wound of debris and neutralise any infectious agents 
that have invaded the tissue. This stage is called inflammation] in a nor¬ 
mal wound inflammation begins a couple of hours after wounding and lasts 
a few weeks. Monocytes metamorphose into macrophages, which complete 
the removal of the pathogens and also secrete some proteins (like Vascular 
Endothelial Growth Eactor, or VEGE). This leads to the next stage: the 
proliferation or reepithelialisation phase. At this point VEGE and other 
substances stimulate the growth and aggregation of the surrounding cells, 
restoring the different tissue functions. The clot is slowly substituted by a 
“temporary skin” called granulation tissue and the interrupted blood and 
lymphatic capillary networks are restored in processes named (blood) an¬ 
giogenesis and lymphangiogenesis^ respectively. After a lengthy period the 
granulation tissue is replaced by normal skin; this happens during the final, 
long phase of remodeling^ which can take up to one or two years. 

Although wound healing has been studied extensively and the main un¬ 
derlying mechanisms are well understood, little is known about how lym¬ 
phangiogenesis takes place. Ear more biological (and mathematical) liter¬ 
ature has been produced about its sibling process, blood angiogenesis; it 
was not until the 1990s that lymphangiogenesis received signiflcant attention 
from researchers [2H1123 EOl- This discrepancy was mainly due to the pre¬ 
vious lack of markers and information on the growth factors involved in the 
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lymphangiogenesis process; such a dearth of biochemical tools impeded a de¬ 
tailed and quantifiable study of lymphatic dynamics 1301111]. For biological 
reviews about lymphangiogenesis see [311111111 and for particularly significant 
biological research papers see mmm- 

Naively, lymphatic vessels may appear “interchangeable” with their blood 
equivalents from a modelling perspective. However, it is stressed that the 
two vasculatures are quite different; for biological papers comparing lym¬ 
phangiogenesis with (blood) angiogenesis see [28l [35l |36]. First of all, the 
capillary structure is completely distinct: while blood vessel walls are rela¬ 
tively thick, surrounded by smooth muscles which pump the blood around 
the body, lymphatic capillaries are made of a single layer of endothelial cells 
known as lymphatic endothelial cells (LEG) [SI]. Moreover, the formation 
of new lymphatic capillaries, or the restoration of preexisting ones, is very 
different from blood angiogenesis. While growing blood capillaries are known 
to sprout from existing interrupted ones, several studies suggest that lym¬ 
phangiogenesis occurs in a different way [231137|. For instance, in |33| it 
is observed that LECs migrate as single cells in the direction of interstitial 
flow and after sufficient numbers have congregated in the wound region, they 
organise into vessels (see Eigurej^. 

We are therefore facing a new process whose mathematical description 
cannot be drawn from any previous model for blood angiogenesis. In the 
following section more details are given about the lymphangiogenesis process 
and a mathematical model is proposed. 

1.3. Outline 

Clearly, successful lymphangiogenesis is an essential element in the wound 
healing process. Yet, as a novel and developing area of attention, its mecha¬ 
nistic basis, relation to other components of wound healing and impairment 
during diabetes remain unclear. With the aim of furthering our understand¬ 
ing of lymphangiogenesis, in this paper a mathematical model is developed 
to describe this process. 

This paper is structured as follows. In Section a brief review of mathe¬ 
matical modelling of wound healing and lymphangiogenesis is given, followed 
by a detailed description of how the terms in our model were chosen. At the 
end of the section, a list of the model parameters and initial conditions is 
provided (a detailed description of parameter estimation can be found in 
. In Section a typical solution of the model is shown and 


Appendix A 
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Figure 1: The photo is taken from [321 Figure 2] and shows lymphatic chan¬ 
nels formation in the mouse tail. Note that at day 10 (a) fluid channeling 
is not observed, but at day 25 (b) discrete channels are present and at day 
60 (c) a hexagonal lymphatic network is nearly complete. Notice also that 
lymph and interstitial fluid flows from left (tail tip) to right (tail base): this is 
in contrast to what happens during blood angiogenesis, which occurs equally 
from both sides of the implanted tissue equivalent. In [32] the authors use a 
new model of skin regeneration consisting of a collagen implant in a mouse 
tail (whose location is indicated by the dashed lines in the picture). The aim 
of the experiment is mainly to characterise the process of lymphatic regen¬ 
eration. Lymph fluid is detected (green in the photograph) and in [32] it is 
shown that LECs follow this fluid. Therefore, this photo can be seen as the 
migration of LECs into the wound. (Bar=lmm) 


compared with real biological datasets. Eurthermore, the diabetic case is in¬ 
troduced and modelled by changing specific parameters. The simulations for 
normal and diabetic cases are then run together and compared with avail¬ 
able data. 


Einally, a steady state analysis (detailed in Appendix B[ ) and 
parameter sensitivity analysis of the model are performed. In Section three 
existing experimental treatments aimed at enhancing lymphangiogenesis are 
presented and simulated; then, potential therapeutic targets are identifled 
based on observations from the parameter sensitivity analysis. Einally, Sec¬ 
tion briefly summarises the main results of the work and possible future 
extensions of the present model are mentioned. 
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2. Modelling 

2.1. Brief review of existing models 

Being of such a complicated nature and evident medical interest, wound 
healing has been the subject of mathematical modelling for decades; see for 
instance [3H1 [39] for the first models in the 1990s, jlD] for a 2002 review and 
m for a 2010 review. A variety of mathematical formalisms have been in¬ 
volved in wound models: from classical PDEs |l2l[33l[3H|, often derived from 
bio-mechanical considerations mms], to stochastic models gS!, to discrete 
models m- Some authors have used moving boundary methods to study 
the movement of the wound edge during healing |1H|, and attempts have 
been made to understand more specific aspects of the healing process, for 
example macrophage dynamics in diabetic wounds [39] and the resolution of 
the inflammatory phase m- One aspect of wound healing that has received 
considerable attention from mathematicians is (blood) angiogenesis (see, for 
instance, [HU IH21 EHl EH EHl EH EH EH])- A comprehensive review of math¬ 
ematical models for vascular network formation can be found in Scianna et 
ahs 2013 review m- Little work has been done with regards to modelling 
lymphangiogenesis and, indeed, m cites only a limited number of math¬ 
ematical works concerning this topic. Of these a representative sample is 
given by isDiini] which deal with tumor lymphangiogenesis and the collagen 
pre-patterning caused by interstitial fluid flow, respectively. More specifi¬ 
cally, m uses the physical theory of rubber materials to develop a model 
explaining the morphology of the lymphatic network in collagen gels, follow¬ 
ing the experimental observations of [32] . This is the only existent model for 
lymphangiogenesis in wound healing known to the authors. A further brief 
review of lymphatic modelling can be found in m, where the phenomenon 
is approached from an engineering perspective. 

In summary, a small number of papers have considered modelling the 
lymphangiogenesis process in the context of tumors; the modelling of lym¬ 
phangiogenesis in wound healing is confined to m, where two fourth order 
PDEs are used to describe the evolution of the collagen volume fraction and 
of the proton concentration in a collagen implant. That work does not ad¬ 
dress the healing process as a whole which is the aim of the present paper. 
Here a simple model (comprised of a system of ODEs) is presented that 
provides an effective description of the main dynamics observed in wound 
healing lymphangiogenesis. 
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2.2. Model development 

In the present model five time-dependent variables are considered: two 
chemical concentrations (TGF-/3 and VEGF) and three cell densities (macro¬ 
phages, LEGs and lymphatic capillaries). Their interactions are described in 
Eigurej^and the formulation of the ODE model is based on the following set 
of processes (the full system is given by the set of equations ®)- The ini¬ 
tial (or pre-wounding) state is altered when latent TGE-/3 is activated (thus 
becoming active TGE-/3, denoted in the sequel by T) by macrophages and 
enzymes released immediately after wounding. This active form of TGE- 
/3 attracts more macrophages (M) to the wound site, through ehemotaxis. 
Macrophages in turn produce VEGE (E), a growth factor that chemoattracts 
and stimulates the proliferation of LEGs (L). Note that LEG growth is also 
inhibited by TGE-^. In the final stage of the process, LEGs cluster in a 
network structure, transdifferentiating into lymphatic capillaries (C). This 
latter process happens spontaneously, although it is enhanced by VEGE. 


In the following the time-dependent variables introduced above are dis¬ 
cussed in detail. In particular, the derivation of the corresponding evolution 
equation is individually presented for each variable. 

TGF-/3 

This chemical is normally stored in an inactive or latent form in the 
body; however, only aetive TGE-/3 plays an important role in wound healing 
lymphangiogenesis, and therefore we will only consider the dynamics of the 
active chemical. Effectively, the active TGE-/3 protein is bound to a molecule 
called Latency Associated Peptide, forming the so-called Small Latent TGE- 
/3 Gomplex. This in turn is linked to another protein called Latent TGE-/5 
Binding Protein, overall forming the Large Latent TGE-/3 Gomplex [63]. 
Hence, the two “peptide shells” must be removed before the organism can 
use the TGE-/3. 

Another feature of this growth factor is that it exists in at least three 
isoforms (TGE-/31, 2 and 3) which play different roles at different stages of 
wound healing |M!- The isoform of primary interest in our application is 
TGE-/31. 

The differential equation describing (active) TGE-/3 concentration has the 
following form: 
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Figure 2: A schematic representation of the model dynamics. The five time- 
dependent variables correspond to the levels of the factors in each of the 
rectangular boxes; the quantity in the circle (latent TGF-/3) is assumed con¬ 
stant. Concerning the arrows, dashed red denotes chemotactic attraction, 
solid green denotes activation/transdifferentiation, dotted blue denotes 
production/enhancement and dash-dotted orange denotes inhibition. 


change in activation by 

TGF-/3 = enzymes and x latent TGF-^ — decay, 
concentration macrophages 

A review of the activation process is presented in [63] where it is reported 
that TGF-/3 can be activated in two ways. The first is enzyme-mediated ae- 
tivation whereby enzymes, mainly plasmin, release the Large Latent TGF-/5 
Complex from the ECM and then the Latency Associated Peptide binds to 
surface receptors. The second form of activation is reeeptor-mediated aetiva- 
tion. Here cells bind the Latency Associated Peptide and later deliver active 
TGF-/3 to their own TGF-/3-Receptors or to the receptors of another cell. 
This behaviour is often observed in activated macrophages [saiM!. 

Thus both enzyme concentration and macrophage density M are influ¬ 
ential in the activation process and thereby appear in the activation term. 
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cipPoe~^p^^^ + Here ap and aM denote the activation rates by en¬ 

zymes and by macrophages, respectively. Tl denotes the (constant) amount 
of available latent TGF-/5 (more details in the next paragraph). In addition, 
the enzyme/plasmin concentration is assumed to decrease exponentially from 
the initial value po, as in ISZl; this reproduces quite well the enzyme dynamics 
in real wounds [68] . 

It is widely accepted that a variety of cells have the potential to secrete 
latent TGF-/3, including platelets, keratinocytes, macrophages, lymphocytes 
and fibroblasts isaEoiins!- Moreover, this latent complex is stored in the 
EGM in order to be constantly available to the surrounding cells [71]. This 
fact is manifested in a constant production rate Tl in our model equation. 
Furthermore, it is well known that macrophages secrete latent TGF-/3 m, 
we assume that this occurs at a constant rate ri. Together these considera¬ 
tions imply that the amount of available latent TGF-/3 in the wound will be 
modelled by Tl + riM. Finally, TGF-^ naturally decays at rate di, so the 
term —diT will be included in the differential equation. Therefore, the full 
equation for active TGF-^ is 
dT 

-jp = [apPo exp{-apTLt) + umM] ■ [L + nM] - diT . ( 1 ) 

Macrophages 

These are a type of white blood cell that removes debris, pathogenic 
microorganisms and cancer cells through phagocytosis. They are produced 
by the differentiation of monocytes and are found in most of the tissues, 
patrolling for potential pathogens. 

Perversely, in addition to enhancing inflammation and stimulating the 
immune system, macrophages can also contribute to decreased immune re¬ 
actions. For this reason they are classified either as Ml (or inflammatory) 
macrophages if they encourage inflammation, or as M2 (or repair) macro¬ 
phages if they decrease inflammation and encourage tissue repair Hence¬ 
forth we restrict attention to inflammatory macrophages, since they are the 
most involved in lymphangiogenesis-related processes. A useful review of 
the multifaceted and versatile role of macrophages in wound healing can be 
found in [73] . 

The following scheme will be considered for macrophage dynamics: 
change in 

macro- _ constant ^ chemotaxis ^ logistic 

phage source by TGF-^^ growth 

density 


removal 
and meta¬ 
morphoses. 
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The various terms appearing in the right-hand side of this equation are dis¬ 
cussed below. 


The number of inflammatory macrophages increases due to their migra¬ 


tion into the wound, in part due to movement of existing inflammatory 
macrophages from the surrounding tissue, as well as by chemotaxis of mono¬ 
cytes up gradients of TGF-^ m, a fraction a of which differentiate into 
inflammatory macrophages The former is modelled by assuming a con¬ 
stant source sm dictated by the non-zero level of inflammatory macrophages 
[76] . and the latter by the term ahi{T) = abiT‘^/{b 2 +T^). Here hi{T) is 
the “chemotactic function”, whose form is discussed in detail in [Appendi^ 
[a} Only a (small) percentage [3 of macrophages undergo mitosis [77|; we 
thus assume the logistic growth term Pr 2 M (1 — M/ki) where r 2 denotes the 
macrophage growth rate and ki the carrying capacity of the wound. No¬ 
tice that here only M appears over the carrying capacity and the other cell 
types L and C are omitted. However, since the logistic term is small overall, 
adding L and C here would just increase the numerical complexity of the 
system without adding any significant contribution to the dynamics of the 
problem. This is reflected in the parameter sensitivity analysis provided later 
in the paper, and simulations (not shown) including all populations showed 
no appreciable difference. Finally, inflammatory macrophages can die, meta¬ 
morphose into repair macrophages or be washed away by the lymph flow. 
This is embodied in the removal term — (^2 + pC)M^ where d 2 denotes the 
constant death rate. Here metamorphosis and removal are considered to be 
linearly proportional to the capillary density C through the coefficient p\ in 
particular, capillary formation is an index of progression through the healing 
process and, to reflect the decreased requirement for inflammatory macro¬ 
phages as wounding proceeds, we assume the metamorphosis/removal rate 
increases with the size of C. Combining these observations one derives the 
macrophage equation 




VEGF 


This is a signal protein whose main function is to induce the forma¬ 
tion of vascular networks by stimulating proliferation, migration and self¬ 
organisation of cells after binding to specific receptors on their surface. There 
are many kinds of VEGF: while VEGF-A and VEGF-B are involved mainly in 
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blood angiogenesis, VEGF-C and VEGF-D are the most important biochem¬ 
ical mediators of lymphangiogenesis via the receptor VEGFR3 (although 
VEGF-G can also stimulate angiogenesis via VEGFR 2 ). For a comprehensive 
description of the growth factors involved in lymphangiogenesis see [7H1E9]. 
Henceforth “VEGF” refers to VEGF-G (and, to a lesser extent, VEGF-D), 
unless otherwise stated. For VEGF we assume the following dynamics: 


change in 
VEGF con¬ 
centration 


constant ^ production by 
source macrophages 


decay 


use by 
LECs. 


Since the normal VEGF level in the skin is nonzero IHOIISI], it is assumed 
there is a constant source of this growth factor from the surrounding 
tissues. VEGF is produced by several cells, but macrophages are considered 
to be one of its main sources in the context of wound healing |82l [83] . It is 
therefore natural to add the production term to the VEGF equation, 

where is the production rate of the chemical by macrophages. On the 
other hand, the VEGF level is reduced by natural decay at constant rate ds, 
taken into account by the term —d^V. In addition VEGF is internalised by 
cells: effectively, LEGs use VEGF to divide and form capillaries [HH[85|; it is 
assumed this process occurs at a constant rate 7 , leading to the term — 7 VL. 
Thus, in this model the equation for VEGF dynamics is 


dV 

dt 


sy + r^M — d^V — 7VL . 


(3) 


LECs 

As discussed above, lymphatic vessel walls are made of (lymphatic) en¬ 
dothelial cells. The equation describing the presence of LEGs in the wound 
consists of the following terms: 



growth, 



transdif¬ 

change in 

upregulated by 

inflow and 

crowding 

ferentia¬ 

LEG 

= VEGF and + 

chemotaxis 

— effect and — 

tion into 

density 

downregulated 

by VEGF 

apoptosis 

capillar- 


by TGF-/3 ies. 


LEG growth is upregulated by VEGF [861 IHIl [85| and downregulated by 
TGF-^ [HHl [89]. The former observation is described mathematically by 
augmenting the normal/basal constant growth rate Ci with V in an increasing 
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saturating manner through the parameters C 2 and C 3 . To account for the 
latter, the growth term is multiplied by a decreasing function of T. Explicitly, 
the whole proliferation term is 



C2 + C3V J 


1 


1 + C4T 


L 


(4) 


where C 4 takes into account the “intensity” of the TGF-/3 inhibition on LEG 
growth. 

It is assumed that LEGs are brought into the wound by lymph flow at a 
constant rate sl and are chemoattracted by VEGE [Ha i. Gonsidering 
a chemotactic function h 2 {V) of the same form as that used for TGE-/3- 
mediated chemotaxis (see Appendix A), these phenomena are captured by 
the terms 

SL + h2iV) = s, + ^^^. 


Since both of these movements originate from the interrupted lymphatic 
vasculature at the edges of the wound, this flow will tend to decrease as the 
lymphatic network is restored. Hence, supposing a linear correlation between 
the term above and the lymphatic regeneration, the former is multiplied by 
the piecewise linear function f{C) defined by 


fiC) 


l-§ iiC <c* 

0 if C > C* 


(5) 


Here C* is a capillary density threshold value above which the lymphatic 
network is functional and uninterrupted and LEG flow stops. Hence the flnal 
term for LEG inflow and chemotaxis is 

LEG growth is limited by over-crowding of the wound space, a fact that is 
taken into account by the negative term —L{M + L + C)/k 2 where k 2 relates 
to the carrying capacity. Finally, individual or small clusters of LEGs migrate 
into the wound and later form multicellular groups that slowly connect to one 
another, organising into vessel structures IH21I311- Here it is assumed that 
when LEGs are sufficiently populous (that is, their density becomes larger 
than a threshold value L*) they self-organise into capillaries at a rate which 
depends linearly on VEGE concentration via the term S 2 V. In particular. 
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endothelial cells tend to form network structures spontaneously (at a constant 
rate (5i, say) but the rate is increased by the presence of VEGF [90]. These 
observations result in the transdifferentiation term — cr(L, C) • {Si + S 2 V)L 
where cr(L, C) is the step function 


a{L,C) 


1 if L + C > L* 
0 if L + C < L* 


(7) 


Note that a depends both on L and C: this is justified by the observation 
that the self-organisation process begins when L reaches the threshold L* 
and then continues as LECs start forming the first capillaries. Therefore the 
complete LEG equation is 

M -\-L-\-C , . . . 

--- L — (j{L^ C) • (5i + S2V)L. 

^2 


m 

( 8 ) 


Lymphatic capillaries 

We assume that the lymphatic capillaries form simply from the self¬ 
organisation of LEGs into a network structure. Thus the capillary formation 
term is just the transdifferentiation term from the LEG equation above and 
the dynamics of C is modelled by 


dt 


a{L,C)‘{5i + 52V)L . 


The full system of equations is therefore given by 


(9) 
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dT 

dt 

dM 

dt 

dV 

dt 

dL 

dt 


dC 

dt 


[ttpPo exp(-apTii) + umM] ■ [Tl + nM] - diT 

activation by latent TGF-^ decay 

enzymes & M$s 

^ +a-^^+Pr2M(^l-^'^- id2j_pCW 

constant ^^ ---^ removal and 

source chemotaxis logistic growth metamorphoses 
by TGF-;d 


sv, + r^M — d^V — 'yVL 

use by LEGs 


constant 

source 


Cl + 


production 
by M$s 

1 / 


decay 

1 


C 2 + ^ 3 ^/ J V 1 ^ C 4 T 


L -\- { Sl -\- 


bsV^ \ 

64 + ^V 


fiC) (10) 


growth upregulated by VEGF inflow and 

and downregulated by TGE-/3 chemotaxis by VEGE 

M -\-L-\-C , . . . 

- - - L — C) • (5i + S2V)L 

k2 ^ ^ -V-^ 

"Z transdifferentiation 

crowding effect into capillaries 

and apoptosis 


(j{L,C)-{5i + 52V)L 

" -V-' 

transdifferentiation of LEGs 


where /(C) and a{L,C) are defined in (j^ and ([^, respectively. 


2.3. Parameters and initial eonditions 
Parameters 

Table [T] gives a full list of parameter values, their units and the sources 
used for their estimation in the normal (non-diabetic) case. It is remarked 
that great care was put into assessing the parameter values, and of the 31 
parameters listed in the table, 25 have been estimated from biological data. 
A detailed description of the estimation of each parameter can be found in 
Appendix A[ 


Initial eonditions 

In the present model, the initial time-point t = 0 corresponds to the 
release of enzymes by platelets within the first hour after wounding [SHlin!. 
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PARAMETER 

VALUE 

UNITS 

SOURCE 

CLp 

2.9 X 10"^ 

mm^pg“^day“^ 

m 

Po 

2.5 X 10® 

pg mm“^ 

no data found 

dM 

0.45 

mm^ cells day 

[HSIISS] 

Tl 

18 

pg mm“^ 

(Ea) 

ri 

3 X 10“® 

pg cells” ^day 

m 

di 

5 X 10^ 

day"® 

m 

sm 

5.42 X 10^ 

cells mm”^day”^ 

(IZSl) 

a 

0.5 

1 

m 

bi 

8 X 10® 

cells pg^(mm^)”^day”^ 

m) 

b2 

8.1 X 10® 

(pg 

miins] 

P 

5 X 10“® 

1 

m 

r2 

1.22 

day"® 

ESI 

ki 

6 X 10® 

mm^cells”^ 

ESI 

d2 

0.2 

day"® 

m 

P 

10“® 

day” ^ cells” ^ 

m 

Sv 

1.9 

cells day”^ 

dSQllHT]) 

rs 

1.9 X 10“® 

pg cells”^day” ^ 

(IH21EH]) 

ds 

11 

day"® 

Eg 

7 

1.4 X 10“® 

mm^ cells ” ^ day ” ^ 

[inni 

Cl 

0.42 


m 

C2 

42 

day 

m 

C3 

4.1 

pg day mm”^ 

m 

C 4 

0.24 

mm^pg“^ 

m 

sl 

5 X 102 

cells day”^ 

no data found 

h 

10^ 

cells pg^(mm^)“^day”^ 

no data found 

bi 

8.1 X 10® 

(pg 

estimated « 62 

G* 

10^ 

cells mm”^ 

m 

k2 

4.71 X 10® 

cells day mm”^ 

m 

L* 

10^ 

cells mm”^ 

m 

<^1 

5 X 10“2 

day"® 

no data found 

<52 

10“® 

mm^pg“^day”^ 

no data found 


Table 1: A list of all the parameters appearing in the model equations (de¬ 
tails of the estimation are provided in Appendix A). Each one is supplied 
with its estimated value, units and source used (when possible) to assess it. 
References in brackets mean that although the parameter was not directly 
estimated from a dataset, its calculated value was compared with the bio¬ 
logical literature; the caption “no data found” signifies that no suitable data 
were found to estimate the parameter. Concerning the VEGF value corre¬ 
sponding to maximum EEC chemotaxis 64, it was assumed that its value is 
similar to its TGE-/3 correspondent this choice was dictated by the lack 
of relevant/applicable biological data, to the authors’ knowledge. 












The initial amounts of active TGF-/3, macrophages and VEGF are taken to 
be their equilibrium values, estimated from experimental data as shown in 
Table It is assumed that there are no endothelial cells or capillaries at 
t = 0. 


INIT.VALUE 

VALUE 

UNITS 

SOURCE 

T(0) 

30 

pg/mm^ 

m 

M(0) 

1875 

cells / mm^ 

m 

y(o) 

0.5 

pg/mm^ 

Eniisi] 

m 

0 

cells/mm^ 

assumption 

C{0) 

0 

cells/mm^ 

assumption 


Table 2: Values of the model variables at t = 0. 


3. Results and analysis 


We now present a typical solution of the system (10) and compare it with 


biological data. The system is solved numerically with the MatLab standard 
ODE solver ode45 with relative tolerance 10“^ and absolute tolerance 10“^ 
over a time interval of 100 days. It is remarked that the present model chiefly 
addresses inflammation and the early proliferation stage of the wound healing 
process. In healthy subjects the inflammatory phase starts a few hours after 
injury and lasts approximately 1 or 2 weeks, but it is prolonged in diabetic 
patients. Moreover, lymphangiogenesis occurs between 25 and 60 days after 
wounding, much later than blood angiogenesis which is observed between 
day 7 and day 17 |2nilSl!- Thus, the equations are expected to realistically 
describe the phenomenon for about the first 100 days post-wounding. 

The TGE-/3 level is expected to display a rapid spike in the first day 
post injury before returning to its equilibrium value [95]. In Eigure the 
simulation output is compared with a biological dataset. Both demonstrate 
the expected initial spike, but in the data a second peak is visible around day 
5, reported also in |9l|. We recall that TGE-^d exists in at least three known 
isoforms: TGE-/31, TGE-/32 and TGE-^3; the biological data set concerns all 
kinds of TGE-^ involved in other wound healing processes, such as collagen 
deposition, which are not modelled here (the time dynamics of the different 
TGE-/3 isoforms can be found in jl()2j ). Nevertheless, the overall predicted 
trend of TGE-^d concentration in the wound matches the biological reality 
fairly well. 
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Figure 3: A comparison of the simulation output for the time course of TGF- 
[3 concentration with data from [95l Figure 2], showing the time course of 
active TGF-/3 generation during wound repair in rats. 


Macrophage levels are observed to reach a peak approximately 5 days 
after injury before returning to their equilibrium level m- Again the model 
prediction is consistent with the biological literature, as can be seen in Fig¬ 
ure 01 



Figure 4: A comparison between model and experimental data for the time 
course of macrophage density. Experimental data taken from |9H Figure 
3c]: note that the time-course comparison here is against the black bars, 
representing macrophage numbers in normal (wild-type) mice. 

VEGF is also reported to reach its maximum concentration 5 days after 
wounding [98]. This is unsurprising given the above macrophage dynamics 
and the fact that macrophages are understood to be primarily responsible 
for the production of the protein. Once again there is a strong correlation 
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between the results of the theoretical model and experimental observations, 
as shown in Figure 


VEGF 




Figure 5: A comparison of the simulation output for the time course of VEGF 
concentration with data from [981 Figure 2], where VEGF was measured in 
rat wound fluid (note that the units on the vertical axis are pg/mL, where 
1000 pg/mL = 1 pg/mm^). 

LEG levels are expected to increase immediately after wounding but only 
later do the LEGs self-organise into capillaries, around day 25 |3l]. This 
is reflected in the simulation shown in Eigure Here LEGs proliferate in 
the wound space until reaching the threshold level L* = 10^ around day 
20. They then start agglomerating in capillary structures, commencing the 
lymphangiogenesis process proper. 

3.1. Modelling the diabetie ease 

In order to simulate the diabetic case, some parameter values are changed 
as described in the following. Unfortunately, it is difficult to obtain precise 
quantitative assessment of the appropriate changes, and therefore the values 
chosen here only have a qualitative significance. 

Several studies report that the TGE-/3 level is significantly lower in dia¬ 
betic wounds compared with controls. This seems to be caused by impaired 
TGE-/3 activation both by platelets and macrophages and by reduced pro¬ 
duction of TGE-/3 by macrophages [1031 11041 I102j . These features of the 
diabetic case are modelled by applying the following modifications to the 
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Figure 6: Simulation output for the time course of LEG density, lymphatic 
capillary density and their sum. Note that the sum density has also been 
plotted the sum, since LEG and capillary cells are difficult to differentiate 
and any cell counts are likely to reflect the total density of these two cell 
types. 


parameters: 

diab _ ^ norm ^ norm diab _ ^ norm ^ norm 

Eurthermore, in diabetic wounds the macrophage density is higher than nor¬ 
mal. In particular, the inflammatory macrophage phenotype persists through 
several days after injury, showing an impaired transition to the repair phe¬ 
notype |l()4j . In addition, macrophage functions (such as phagocytosis and 
migration) are impaired in the diabetic case |105lll()6| . These differences from 
the normal case are reflected in the following choice of parameter changes: 


^diab _ 

0.8 > a "'”’’”, 

Jjdiab 

^diOib _ 

1 Jnorm ^ ^norm 

2^2 ^ ^2 5 

r^diab 

'3 


3 imorm ^ imorm 

4^1 ^ 5 


l^norm ^ ^norm 
2'3 ^ '3 


Finally, it is reported that endothelial cell proliferation is markedly re- 


19 







duced in diabetic wounds when compared with the normal case [TTITlITMrmn] 
(a detailed discussion of endothelial dysfunction in diabetes can be found in 
mni). This phenomenon is reflected in the model by reducing the basal 
proliferation rate of endothelial cells: 

diab _ ^ norm ^ norm 

Cl ~ ^ 

3.2. Comparison of results in the normal and diabetie eases 

Figures [7jjT^ show numerical simulations of the model comparing the time- 
course of the flve model variables in the normal (blue solid line) and diabetic 
(red dashed line) case. 

TGF-f3 

The level of TGF-/3 recorded in diabetic wounds is lower than in the 
normal case, at least in the first 20 days after injury [inSlIinillinS]. Model 
simulations are consistent with this (Figure [^. 


Active TGF-p 





Figure 7: Time course of TGF-/3 concentration in normal and diabetic 
wounds as predicted by the model and as found empirically in |102l Fig¬ 
ures 3G and 3H], where the authors study molecular dynamics during oral 
wound healing in normal (blue lines) and diabetic (red lines) mice. 
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Macrophage 

Experiments show that the density of macrophages in diabetic wounds is 
higher than in the normal case and they persist for longer in the wound site 
|104[ [731 Il06| . Model simulations match these observations (Figure [^. 


Macrophages 




Figure 8: Time course of macrophage density in the normal and diabetic 
cases. We compare the model prediction with jlQ4[ Figure IB]. In the ex¬ 
perimental results, the relative height of shaded to solid white bars indicates 
the relative macrophage density in diabetic/non-diabetic wounds in mice 
(assessed via Ly6C expression, a marker for the macrophage lineage). 


VEGF 

The VEGF level during wound healing is lower in diabetic patients una 
llQ4j . In fact, as described below, a key target for the design of new therapies 
has been increasing VEGF levels. The simulation output and a biological 
dataset are compared in Figure 

LECs & Capillaries 

In diabetic patients, lymphatic capillary formation is delayed and insuf¬ 
ficient uni Ea I21!- The model simulations are consistent with this (Figure 

T§. 

3.3. Analysis of the model 

In this section the steady states of the system are identified and a sensi¬ 
tivity analysis of the model parameters is performed. 
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VEGF 



□ C 



Time (Post wounding) 


Figure 9: Time course of VEGF in normal and diabetic case wounds. Com¬ 
pare the simulation with the data reported in jlOSi Figure 2G], represent¬ 
ing VEGF presence in control (white bars) and diabetic (black bars) rats. 
In fact, investigates the connection between a defect in resolution of 

inflammation and the impairment of TGF-/3 signaling, resulting in delayed 
wound healing in diabetic female rats. The abbreviations in the legend stand 
for: C, control; D, diabetic; +E2, diabetic-treated with estrogen; TNFRl, 
diabetic treated with the TNF receptor antagonist PEG-sTNF-RI; VEGF, 
vascular endothelial growth factor. 


3.3.1. Steady States 

For the parameter set studied, at t = 0 there are no LECs in the wound, 
but subsequently they increase towards a positive value of approximately 
2 X 10^ cells/mm^. However, when they reach the “threshold” density L* = 
10^ cells/mm^, the system steady states change and L starts to decrease 
towards zero. In the meantime, lymphatic capillaries start forming; their 
final value will depend on the dynamics of the system, but in any case it will 
be bigger than C* = 10^ cells/mm^. This “switch” in the steady state values 
is due to the presence of two piecewise defined functions {a and /) in the 
system. On the other hand, there is always one stable steady state for M 
which also defines one for T: 

1875 cells/mm^ and ^ 30 pg/mm^ . 

di 

Note that the M steady state (and thus also that for T) is unique for param¬ 
eters with biologically relevant values. Eor the V-equilibrium, the following 
expression is found: 

yeq ^ Sy + 

4 + 7 ^"'' 
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Figure 10: Time course of LEG and capillary density in normal and diabetic 
cases, compared with [231 Figure 2b]. The study presented in [23] investigates 
the role of wound-associated lymphatic vessels in corneal inflammation and 
in a skin wound model of wild-type and diabetic mice. The flgure shows a 
quantification of lymphangiogenesis in the corneal suture model assay in the 
wild-type (db/+) and diabetic (db/db) cases. No suitable data were found 
in the skin wound model. 


Note that depends on L; therefore V will tend to different values accord¬ 
ing to the current L-steady state; for = 0 it is — 0.5 pg/mm^. Details 
about how these steady states were determined can be found in [Appendix 


The stability of the steady states is determined numerically. The stability 
of is deduced from the shape of the numerically-plotted M-nullcline, and 
the stability of the other steady states can be inferred from the simulations 
of the full system. See for instance the simulation shown in Figure pT| where 
the model is run over a time interval of 250 days: here it is evident that 
all the variables tend to stay at a stable value after about 100 days post- 
wounding. Since some of the parameters were modifled to simulate diabetes- 
related conditions, the steady states for the diabetic case are different than 
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the corresponding “normal” ones. In particular, TGF-/3, VEGF and capillary 
equilibrium values are lower in the diabetic case, while the macrophage level 
is higher than in the normal case. LEGs go to zero in both cases. 

Although this analysis does not give very profound insights into the un¬ 
derstanding of lymphangiogenesis, it provides some extra information about 
the dynamics of the system. More specifically, it shows that for realistic pa¬ 
rameter values the system has only one steady state, which is in agreement 
with experimental observations. 


3.3.2. Parameter Sensitivity Analysis 

Here a numerical parameter sensitivity analysis of the model is presented 
which plays two important roles. On the one hand, it demonstrates which 
parameters are most significant in the model, and thereby provides a deeper 
understanding of the dynamics of the system. On the other hand, it consti¬ 
tutes the first step towards the design of new therapeutic approaches. 

To estimate the dependence of the model on a given parameter p, a quan¬ 
tification of the affect of a change in p on the (final) capillary density C at day 
100 is calculated. To begin, p is increased by 10% and thereafter the system is 
solved over the time interval [0,100]. The final value of the capillary density 
thus obtained, denoted is then compared with the reference value 

of the corresponding density in the original system. The percentage 
change is defined by 


percentage change in C = 100 x 


^+ 10 % _ (jref 


C^ef 


( 11 ) 


The same procedure is then repeated, this time substituting the parameter 
p with its value deereased by 10% and the corresponding change is 

calculated. The results are summarised in Eigure 

It is notable that perturbing any parameter does not result in a percentage 
change in final capillary density of more than 15%. Thus, the model is quite 
robust in terms of dependence on the parameters. Percentage changes over 
5% are observed only for eight parameters. Of these, one needs to deerease 
dMi Tl^ sm or C 4 to observe an increase in the final capillary density; while 
a similar enhancement is obtained by inereasing di, ^ 2 , or k 2 . 
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Figure 11: Simulation of the model in both normal (solid blue) and diabetic 
(dashed red) cases over a time period of 250 days. 
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Figure 12: Percentage change in the final capillary density C(IOO) when every 
parameter is increased/decreased by 10%. 


4. Therapies 

J^.l. Existing experimental treatments 

Although there is at present no approved therapy for enhancement of 
lymphangiogenesis (in wound healing or in any other context), many studies 
and experiments have been published exploring potential treatments. In the 
following three such experiments are reported and then simulated. 

Administration of TGF-P Reeeptor-Inhibitor 

This substance binds to the TGF-/3 receptors on the surface of surround¬ 
ing cells, thus making them “insensitive” to TGF-/5 molecules and their effect, 
mu reports a study of the effect of TGF-/3 on lymphangiogenesis in which 
human LEGs are cultured and quantified after treatments with TGF-/31 or 
T/3R-I inhibitor to assess cell growth, cord formation and cell migration. It is 
observed that TGF-^dl treatment decreases cord formation, while T/5R-I in¬ 
hibitor treatment increases it. These results are consistent with those found 


26 































































cells/miT? pg/mm' 


in m, where it is reported that a higher level of TGF-/31 is associated 
with delayed recruitment and decreased proliferation of LECs during wound 
repair. 

To simulate the treatment with T^dR-I inhibitor, the cell migration assay 
is considered. Here, the inhibitor was added at 3 fiM = 817 pg/mm^ (the 
molecular weight is 272). Since this is significantly bigger than the concentra¬ 
tion of TGF-/5 in our model and in normal skin (in both cases the maximum 
level is 300 pg/mm^), this treatment is simulated by setting the parameters a 
and C 4 equal to zero (that is, TGE-^d molecules have no effect on cells because 
their receptors are “occupied” by the inhibitor). The effect of this “virtual 
treatment” are shown in Figure and match well with the described TGF- 
/3 inhibitor experiment: LEG and capillary densities are markedly increased 
by the treatment. 
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Figure 13: Time courses of T, M, F, L, C and L + C in a simulation of the 
T/5R-I inhibitor experiment described in m- 


Macrophage-based treatments 

Another therapeutic approach is to add macrophages to the wound, so 
that they secrete VEGF and other substances that are known to induce lym- 
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phangiogenesis. In m an “opposite” experiment is implemented: here a 
systemic depletion of macrophages is reported to markedly reduce lymphan- 
giogenesis. This is in accordance with m, in which it is observed that the 
induction of macrophage apoptosis inhibits IL-l^-induced lymphangiogene- 
sis. One hypothesis suggested to explain such results is that because of the 
reduced level of macrophages, less VEGF is produced and this impairs LEG 
proliferation and capillary formation. 

We simulated the increase in macrophage apoptosis by taking a bigger (for 
instance, the double) value of ^2 in the system. The output of the model in 
which ^2 is doubled (both in normal and diabetic cases) is reported in Figure 
14} In this case, the output is in contrast with what is described in the 


biological studies: although fewer macrophages and consequently less VEGF 
are present, more LEGs and capillaries form after the simulated treatment. 


Active TGF-p 
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Figure 14: Time courses of T, M, V, L, C and L + C in a simulation of the 
macrophage-depletion experiment described in |113| . 


This result could be explained by the fact that, in the model, a reduction 
in macrophage density implies a reduction in TGF-/5 level, so that the inhi¬ 
bition of LEG proliferation is smaller and hence there are more endothelial 
cells to form the capillaries. In fact, in the previous section it was found that 
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the system is much more sensitive to C 4 than to C 2 , C 3 or 82 . It is then natural 
to consider the effect of fixing T — 30 pg/mm^ in the LEG growth term Q; 
this level of T corresponds to the TGE-^d equilibrium. The simulation output 
in this case is shown in Figure 
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Figure 15: Time courses of T, M, F, L, C and L + C in a simulation of 
the macrophage-depletion experiment described in |113j . where the T in the 
LEG growth inhibition term is substituted by = 30 pg/mm^. 


With T fixed, the difference between the treated and untreated cases is 
very small, but still an increase in capillary formation is observed in spite of 
the lower VEGF level. This may be due to the fact that, with fewer macro¬ 
phages, the crowding term in the LEG equation is smaller, which facilitates 
the growth and accumulation of endothelial cells. In fact, if the M in the 
crowding term is fixed at its equilibrium value of 1,875 cells/mm^, there is no 
difference at all between treated and untreated cases. Note that this result 
could have been foreseen from the parameter sensitivity analysis (Figure 
which predicted that a 10% increase in d 2 induces a 5 to 10% increase in 
final capillary density. 
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VEGF supply 

A third documented approach to enhance lymphangiogenesis consists of 
supplying VEGF to the wound, since this protein promotes both LEG growth 
and the ability of LEGs to form a network-like structure. For instance, in 
m a wound healing assessment is done in normal and diabetic mice after 
a VEGF-treatment. More precisely, two different types of VEGF were stud¬ 
ied: VEGF-Ai 64 and VEGF-Etvzt- The authors observed that the treatment 
with VEGF-Ai 64 increased macrophage numbers and the extent of lymphan¬ 
giogenesis in both wild-type and diabetic cases, while VEGF-Etvzt does not 
induce any significant change. 

In order to reproduce the experiment in silico^ the amount of supplied 
VEGF is estimated in [Appendix A.4[ Then, a 10 days constant VEGF 


supply of 1.8 X 10^ pg/mm^ is introduced in the model system. The output 
is reported in Figure \TE\ 
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Figure 16: Time courses of TGF-^, macrophage, VEGF, LEG and capillary 
densities in a simulation of the 10-days VEGF-supply experiment described 
in Here the original model is altered by adding a constant VEGF 

supply of 1.8 X 10^ pg/mm^ for 0 < t < 10. 


There is apparently no difference between capillary formations of treated 
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and untreated cases. Moreover this result is relatively insensitive to the 
amount of VEGF supplied. What if the same treatment is applied for 30 
days instead of 10? A simulation of this is shown in Figure [T^ 
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Figure 17: Time courses of M, L, C and L + C in a simulation of the 30-days 
VFGF-supply experiment. 


There is now a clear difference in the treated cases (especially the diabetic 
one) showing a lower level of macrophages and an earlier onset of capillary 
formation, even if the final capillary density is similar to that in the untreated 
case. 

Observation 

There is an important feature common to all three modelled therapies: in 
order to stimulate lymphatic capillary formation one cannot consider TGF- 
/3 or VFGF levels individually. A precise balance of TGF-/3 and VFGF is 
necessary for successful lymphangiogenesis. This mutual equilibrium may 
be reached in vivo by the production of either of these growth factors by 
other cell types not considered in the present work. In particular, this sug¬ 
gests the model does not take into account certain elements of the process. 
Nevertheless, the model does effectively describe both normal and diabetic 
lymphangiogenesis in wound healing which suggests that the variables con¬ 
sidered here are the most relevant. This indicates that potential therapies 
should focus on these aspects of the regeneration process. 

4^2. Novel therapeutie approaehes 

As mentioned above, parameter sensitivity analysis proves useful in de¬ 
signing novel therapeutic approaches. Among the “sensitive” parameters 
only um, ^2 and Ci vary between the normal and diabetic cases (note that 
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the increasing d 2 -case was discussed above in the macrophage-based treat¬ 
ment). Thus, at least theoretically, aM, <^2 and Ci are the natural targets 
for a therapeutic strategy aiming to increase the final lymphatic capillary 
density. The feasibility of each suggested parameter change is now explored. 

Decreasing au 

Decreasing au means lowering the macrophage-mediated activation of 
TGF-/5. First of all, note that the increase in final capillary density due 
to a decrease in au is explained by the fact that less active TGF-/3 implies 
less TGF-^-inhibition of LEG growth and hence a larger LEG growth term. 
Eor the practical implementation of this change, it is recalled that receptor- 
mediated TGE-/3 activation consists of the binding of Latency Associated 
Peptide to the cell surface through receptors such as TSP-1/GD36, M6PR 
and multiple aV-containing integrins |63]. Hence, a decrease of aM might be 
obtainable by blocking these receptors. 

Increasing Ci 

Increasing ci would be achieved by increasing the LEG growth rate. Sev¬ 
eral possible implementations of this are found in the literature. Recombi¬ 
nant human IL-8 induces (human umbilical vein and dermal microvascular) 
endothelial cell proliferation and capillary tube organization m- DNA 
dependent protein kinase (DNA-PK) is well known for its importance in re¬ 
pairing DNA double strand breaks; in m\ it is observed that DNA-PKcs 
suppression induces basal endothelial cell proliferation. In HH! it is reported 
that poly dopamine-modified surfaces were beneficial to the proliferation of 
endothelial cells. Einally, non-thermal dielectric barrier discharge plasma is 
being developed for a wide range of medical applications, including wound 
healing; in particular, m reports that endothelial cells treated with plasma 
for 30s demonstrated twice as much proliferation as untreated cells, five days 
after plasma treatment. 

Other parameters 

To increase the final capillary density, one could also think about targeting 
other parameters to which the system is sensitive. In particular: 

• Decreasing Tl means reducing available (latent) TGE-^ and hence re¬ 
ducing the TGE-^-inhibition over LEGs. Suppression of TGE-/3 by 
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antibodies has been proposed as a possible therapy to reduce scar for¬ 
mation jl 2 ()l 11211 I 122 j . Thus many studies of TGF-^ antibodies are 
available. 

• Decreasing C 4 involves reducing the (inhibitory) effect of TGF-/3 on 
LEGs, which can be achieved by blocking specific TGF-/5-receptors on 
the endothelial cell surface. Now, TGF-^ signalling is very well studied 
[T2S!; in particular, it is known that TGF-/3 family proteins act through 
two type 11 and two type 1 receptors and that ALK -1 antagonizes the 
activities of the canonical TGF-/3 type 1 receptor, T/3R1/ALK-5, in the 
control of endothelial function HM]. Moreover, a few studies have been 
published which deal with blocking of TGF-/3 receptors in the specific 
case of endothelial cells fmfTMfTTF] . 

• Ghanges in the other “sensitive” parameters do not appear feasible. 
Increasing di would mean increasing the TGF-^d decay rate; decreasing 
sm would mean reducing the constant source of macrophages; increas¬ 
ing k 2 requires an increased “carrying capacity” for the wound. We are 
not aware of practical approaches that could cause these changes. 

• Finally, among the parameters that, when changed by 10% of their 
value, induce a change in final capillary density between 2 and 5% 
(that is, a bit less than those analysed above), only bi merits discussion. 
Reducing hi corresponds to reducing macrophage chemotaxis towards 
TGF-^, which might be achievable by blocking specific receptors on 
the macrophage surface. 

5. Conclusions 

Our model procures new insights into the mechanisms behind lymphan- 
giogenesis in wound healing. The major contributors to the process have been 
identified (TGF-/3, macrophages, VEGF and LEGs); the self-organisation hy¬ 
pothesis for the lymphatic network formation described in [3211311 has been 
confirmed and the importance of the balance between different factors has 
been highlighted. Moreover, the present work suggests novel therapeutic ap¬ 
proaches to enhance the lymphangiogenic process in impaired wound healing. 
In addition, nearly all of the relevant parameters have been estimated from 
biological data and therefore this work provides fairly reliable numerical val¬ 
ues for the parameters encountered. However, any parameter estimation is 
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limited by, for example, the specific experimental method used or discrepan¬ 
cies between the system considered here and that studied in a given reference. 
The results should therefore be viewed with care. In particular, the numeri¬ 
cal values pertaining to the aforementioned balance between the TGF-^ and 
VEGF may be shifted under an alteration of the parameter set. 

This paper is intended as a first step in studying wound healing lymphan- 
giogenesis through mathematical modelling. Further work should include a 
spatial variable (and thus involve PDFs rather than ODEs) in order to take 
into account the important role of lymph flow in lymphatic capillary network 
formation. Introducing a spatial variable would also enable a fuller descrip¬ 
tion of chemotaxis. A PDE model would also be able to reflect further differ¬ 
ences between angiogenesis and lymphangiogenesis. In particular, contrary 
to blood angiogenesis, lymphangiogenesis is unidirectional: as opposed to 
sprouting from both sides of the wound, LEGs appear to predominantly mi¬ 
grate downstream to the wound space in the direction of the interstitial flow 
m The model could also be extended to include other aspects of wound 
healing, such as blood angiogenesis: implemented effectively, this would give 
a more detailed overview of the different mechanisms and the time-scales 
involved in the various processes. 
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Appendix A. Parameter estimation 

Appendix A.l. Equilibrium values and standard sizes 
TGF-f3 equilibrium 

The equilibrium value of active TGE-/3 is about 30 pg/mm^ [95l Figure 

2 ]. 

Maerophage equilibrium 

The macrophage steady state can be estimated from m Figure 1], which 
plots typical macrophage density in the skin. This shows that there is an 
average of about 15 macrophages per O.lmm^ field. Assuming a visual depth 
of 80 /im, the macrophage density becomes 15 cells/(0.1mm^ x 0.08mm) = 
1875 cells/mm^. 
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VEGF equilibrium 

The VEGF equilibrium concentration is estimated to be 0.5 pg/mm^ from 
[801 Figure 1] and jSU Figure 2]. 

Final LFC and Capillary density 

In [31] we find that “it was not until day 60^ when functional and contin¬ 
uous lymphatic capillaries appeared normal” and “at day 60 the regenerated 
region had a complete lymphatic vasculature, the morphology of which ap¬ 
peared similar to that of native vessels”. Hence, if we assume that a capillary 
network that can be considered “final” appears at day 60, we will take 
(or to be the number of LECs present at this time, and since in the nor¬ 
mal case all the lymphatic endothelial cells will become part of the capillary 
network, we will further assume that — 0. In |3H Figure 2E] we see that 
at that time there are about 80 cells. This value corresponds to a 12 jam thin 
section. In addition, through PI Figure 2D] we can calculate the observed 
wound area, which is about 5.6 x 10^/xm^. In this way we get a volume of 
0.0067 mm^ with 80 cells, which corresponds to = 1.2 x 10^ cells/mm^. 

FC size and weight 

[128] reports the cross-sectional area of an EC as lO/xm x 100 jam. In [129] 
the thickness of these cells is given as 0.5/im. Hence we can assume a cell 
volume of approximately 500/im^ = 5 x 10“^ mm^. Moreover, if we assume 
the density of the cells to be 1 g/mL (the same as that of water), we have 
that a cell weighs about 5 x 10“^° g = 500 pg. 

VEGF moleeular weight 

The molecular weight of VEGF is 40kDa = 66.4 x 10“^ pg/mol [130] . 

Appendix A.2. TGF-^ equation 
Enzyme-mediated aetivation rate a^ 

For ttp, it seems reasonable to take the rate at which LAP binds to the 
receptors. Now, in m we find an estimate for the binding rate to be about 
1.7 X 10^ M“^s“^ [9ll Tables I and IV]. Considering a protein weight of 
approximately 50 kDa (found in the same article) and converting the units 
we find: Up ^ 2 x 10“^ mm^pg“^min“^ = 2.9 x 10“^ mm^pg“^day“^. 
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Receptor-mediated activation rate om 

[66] reports that (“activated”) macrophages plated at 2 x 10^ cells/well 
in a 24-well tissue culture dish (that is, about 10^ cells/mm^) activated ap¬ 
proximately 8% of the total TGF-^ secreted after 22 hours. This means that 
1 cell/mm^ activated about 0.087% of TGF-^ per day. On the other hand, 
in [65] we find that macrophages cultured at 1 x 10^ cells/dish in 100-mm 
plates (that is, about 10 cells/mm^) activated 12.2% of total TGF-/31 after 
about 36 hours. So 1 cell/mm^ activated approximately 0.8% of the TGF-/3 
per day. Hence, for om we will take a value between 0.087 and 0.8, say the 
average 0.45. 

TGF-f3 production rate ri 

In [70] it is reported that 10^ macrophages produced about 30 pg of TGF- 
/3 after 24 hours. So one single macrophage produced 30 x 10“^ pg/day of 
TGF-^, and we can then take ri = 3 x 10“^ pg • cells“^ • day“^. 

TGF-P decay rate di 

In |93| it is stated that “free TGF-/3 has a half life of about 2 min”. We will 
therefore take the decay rate di of active TGF-^ to be di — In 2/(2 min) = 
0.35 min“^ = 500 day“^. 

Gonstant amount of latent TGF-^ Tl 

At the equilibrium, the TGF-/3 equation becomes: [0 + auM^^] • \Tl + 
riM^^] — diT^^ — 0 . Substituting the values of the parameters um, di and 
of the equilibrium values and found before, we get an equation for 
Tl. Solving it, we find Tl — 18.0916 pg/mm^. To compare this value with a 
“real” one, we consider m- Here, taking an average of 6 pg of TGF-/31 per 
mg of skin (from [92( Figure 3]) and assuming a skin density of 1 g/mL (as 
for water), we have a concentration of latent TGF-^d of 6 pg/mm^, which is 
of the same order as our previous estimate. 

Appendix A.3. Macrophage equation 

Fraction of monocytes migrating into the wound that differentiate into macro¬ 
phages a 

We follow |19| and take a to be equal to 0.5 in normal wound healing 
(reflecting the fact that in this case the number of inflammatory macrophages 
is the same of the repair ones), and a = 0.8 in a diabetic wound (since there 
are more inflammatory than repair macrophages this time). However, it must 
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be noted that m comments that “there is currently no quantitative data on 
which the value of a for diabetic wounds can be based”. 


Migration of monocytes to the wound in response to TGF-f3: shape of hi{T) 
and value of 62 

In [71] the authors study the migration of monocytes taken from healthy 
volunteers and observe that the motility of the cells depends on the dose 
of TGF-/5 to which they are exposed. Moreover, even if the response varied 
with individual donors, they see that “the optimal chemotactic concentration 
for TGF-/5 fell within the range 0.1-1.0 pg/mL”. These findings are shown 
in Figure A. 18, which shows the dataset reported in m In light of these 
observations, it is reasonable to take the chemotactic function hi{T) to be 


h{T) = 


biT 


(62 + 


where we will take m — 2 from visual comparison with data. 

To determine 62 we look for the maximum of hi{T). This is located at the 
point Tmax where h'^{Tjnax) — 0 , which is T^ax — [HS] reports that the 

normal level of TGF-^d in the skin is about 30 pg/mm^, and that this amount 
increases up to 300 pg/mm^ during wound healing. From this, and from the 
observation that macrophage levels also reach a peak soon after this TGF-/3 
peak |94|, we deduce that the maximum monocyte/macrophage migration 
occurs when the level of TGF-/3 in the skin is about 300 pg/mm^. Thus we 
take 62 = 300^^; for m = 2 this gives 62 = 300^ = 8.1 x 10 ^ pg^ (mm”^)"^. 
This seems to be in contrast with the data reported in [TH Figure 1]: here 
Tjy^ax is around 0.5 pg/mL = 0.5 x 10“^ pg/mm^. However, this value does 
not seem realistic and should be considered carefully. In particular, we recall 
that chemotaxis occurs through gradients of a chemical, and considering 
this directed movement to depend only on the absolute concentration of the 
chemoattractant is a simplification. In fact, in the experiment described 
in m the diluted chemotactic stimuli were placed in the bottom wells of 
microchamber plates that were separated from the upper wells by a filter 
with 5.0 fim pores; then monocytes were put in the plates with human TGF- 
[3 diluted at different concentrations. Ghemotactic activity was defined as 
the mean number of monocytes that migrated through the pores. Therefore, 
we consider the value of T^ax fo be a “conventional” one and we keep the 
estimate for 62 found above. 
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Figure A.18: Experimental data from m Figure 1] reporting a quantification 
of monocyte chemotaxis for different concentrations of TGE-^d. Chemotactic 
activity is defined as the mean number of monocytes that migrated through 
the 5-/xm pores in three standard fields for each of triplicate filters. 


Percentage of monocytes/macrophages undergoing mitosis / 

In [771 Figure 4] we see that only about 0.5% of monocytes in vitro show 
mitotic activity. Therefore we take / = 0.005. 


Macrophage growth rate r 2 and carrying capacity ki 

To obtain estimates for these parameters, we first consider just the logistic 
part of the macrophage equation dM/dt — r 2 M(l — Mjk/)^ whose solution 
is 


M{t) 


A;i + Mo(e"2t_i) ' 


(A.l) 


In [96] murine macrophage-like cell growth is measured under different con¬ 
ditions. We then fit the data contained in [961 Figure 1] to the curve 
) using the MatLab function nlinfit. Moreover, [96] specifies that 
“6 X 10^ cells were cultured in 100-mm tissue culture plates in 10 mL of 
the (above) medium”. Then, taking Mq = 600 cells/mm^, we get the esti¬ 
mates f2 = 1.22, fci = 6 X 10^, with 95% confidence intervals (0.82,1.62) and 
(4.40 X 10^,7.61 X 10^) respectively (these are calculated using the MatLab 
function nlparci). 
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Macrophage constant removal rate ^2 

m presents a mathematical model for keloid and hypertrophic scarring. 
Here we find that macrophages are known to exist in a wound for several 
days after the initial migration; so based on this we assume a decay rate for 
macrophage cells to be of the order of ^2 ~ 0.2 day“^. 

Macrophage capillary-dependent removal rate p 

We assume that the term pC becomes of the same order of ^2 when 
capillaries reach their “final” density, which we have estimated as 1.2 x 10^. 
Since we estimate d 2 — 0.2, we take p — 10“^. 

Migration of monocytes to the wound in response to TGF-^: bi 

To find bi we notice that this parameter determines the maximum level 
of macrophages M during healing. To find this value, we refer to m Here 
the authors investigate the role of TGF-/5 activation in wound repair, and 
assess several components of wound healing (including inflammatory cell in¬ 
filtration) over a period of 28 days. [9lj reports that at day 5 a maximum 
of about 70 macrophages/field (400x) are observed. Assuming a diameter 
field of view of 0.4 mm and a depth of field of 80 /xm, 1 field corresponds 
to (0.2)^7r mm^ x 80 x 10“^mm 10“^mm^. Then 70 cells/field 7000 

cells/mm^. Numerical experimentation shows that reproduction of this re¬ 
sult requires 6i 8 x 10^ cells mm pg^/day. 

Macrophage constant source sm 

At the steady state, the M-equation becomes 

SM + + /5^2M (^1 - - d2M - pCM = 0 . 

Substituting the equilibrium values and the parameters found above, we get 

Sm = 586 - X 5.5 x 10“^ . 

Above, we chose bi — 8x10^ cells mm pg^/day. Therefore sm — 542 cells/day. 

Appendix VEGF equation 

VEGF production by macrophages r^ 

From IH21 Figure IB] we have that human macrophages plated at 10^ cells/ml 
lO^cells/mm^ produced 214 pg/ml = 214 x 10“^pg/mm^ of VEGF after 24 
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hours of culture. Then rs 8.9 x 10“^pg-cells“^-h”^ = 2.1x10“^ pg-cells“^- 
day“^. In fact, this value is a bit smaller than the one we use in the model, 
which we take to be 9 x 2.1 x 10“^ ~ 1.9 x 10“^ pg/cell/day. This “adjust¬ 
ment” is done considering the data shown in [98l Figure 2], which reports 
that the VEGF peak (occurring at day 5) corresponds to a level of about 
1000 pg/mL = 1 pg/mm^. Since in our model we assume that the VEGF 
peak is due (mainly) to the macrophages, it seems reasonable to adjust the 
parameter rs to meet this observation. 

Other estimates were obtained from IHS1II3I] and one “equivalent” param¬ 
eter was found in the modelling paper jl32j . Although the numerical values 
are different, they are all between 10“^ and 10“^ pg/cell/day. This variety is 
not surprising because different cell types produce VEGF at different rates 
(as clearly shown in jl31j ). Since in the context of wound healing lymphan- 
giogenesis we are mainly concerned with macrophages, we focus more on the 
values for these cells [821 [83] . 


VEGF decay rate 

The half-life for VEGF 165 (the most common and biologically active 
VEGF protein) at room temperature is 90 minutes [99]. It follows that 
ds = 11 day“^. To compare this value with those used in other modelling 
articles, we see that in the VEGF decay rate is taken to be jiy — 

0.456 h“^ = 10.9 day“\ while in [58] the VEGF natural decay/neutralization 
rate is /ic = 0.65 h“^ = 15.6 day“^. 


VEGF consumption/internalization by EGs 7 

In [lOOj the VEGF internalization by a cell is described in the way schema¬ 
tised in Figure [AT^ In that figure and the following text, V stays for VEGF, 
R for receptor and VR for the ligand-receptor complex. 

We will not distinguish between the two kinds of receptors R1 and i? 2 , so we 
will just take the average of the above parameters. That is, we will take 

fcon = 2.5 X 10®M“^S-1, = 2.5 x IQ-V^, . . 

= 10 = 2.8 X 10 ^. 


According to this reaction scheme, the corresponding equation for V (for a 
single cell) is 


d^ 

dt 


= —/Con • E • i? + /Coff ‘ (VR) . 


(A.3) 
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PARAMETER VALUES 



kony^Ri = 3.8 X 
^on,y,R2 — 1-2 X lO^M ^ 

^off,v,Ri = 95 X 10-^s-i 

^off,y,R2 ~ "^^9 X 10 ^ 

^int,Ri “ ^int,R2 = 1 x 10 
^int,VRi ^ ^int,VR 2 = 28 x 10 


Figure A. 19: A schematic representation of the dynamics considered in 
V stands for VEGF, R for Receptor and VR for the ligand-receptor complex. 


Here the dimensions of R are moles per unit volume. Similarly the equation 
for the ligand-receptor complex (VR) (per cell) is 

= +kon -V-R-k^s- (VR) - ' (VR) ■ 

So, at equilibrium: 


/Con ‘ V ‘ R — /CqIJ • (VR)eq — ‘ R)eq — 0 

^on ‘ R 


{VR)eq = 


^Off + hntyR 


V 


(A.4) 


To determine the value of i?, we note that human endothelial cells (cultured 
in vitro) display 1,800 VEGFRl/cell and 5,800 VEGER2/cell jl34j . giving a 
total of 7.6 X 10^ receptors/cell. Assuming a cell volume of 500 /xm^ we take 
R = 2.5 X 10“^^ mol/mm^. Substituting this value for R and replacing the 
parameters with the values given in (A.2), the equation (A.4) gives (VR) 


eq 


118 • E = VReq • V. Substituting this value for (VR) in (A.3) we have 


dV 

dt 


~ ’ di ~ ^off ‘ di^q^ ‘V ‘ L 

=2.9x103 day“^ 


Then, again assuming a cell volume of 500/im^, we have that our parameter 
(per cell) is 

7 = 5 X 10“^ mm^ x 2.9 x 10^ day“^cells“^ = 1.4 x 10“^ mm^ • day“^ • cells“^. 
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VEGF constant source Sy 

At equilibrium, the ^/-equation becomes 


+ r^M — d^V — ^VL = 0 


Substituting the equilibrium values and the parameters gives = 1.9 cells/day. 
VEGF supply in 1115^ 

In jll5j control medium with or without growth factors was injected into 
wounded mice at a dose of 2 /ig/wound/10 days. Since the experimental 
wounds are full-thickness with a 5-mm diameter, and assuming a skin thick¬ 
ness of 0.56 mm unn], this amount corresponds approximately to 1.8 x 10^ 
pg/mm^/day. Now, since the experiment is performed in vivo^ it is reasonable 
to assume that the vast majority of the VEGF is washed away and dispersed 
by body fluids. In the absence of quantitative data, we assume that only 1% 
of the added VEGF is therapeutically active in the wound, giving a delivery 
of 1.8 X 10^ pg/mm^/day. 


Appendix A.5. LEGs equation 

We consider first the parameter Ci, then A; 2 , and then the remaining pa¬ 
rameters. For Cl and k 2 we focus on the “logistic” part of the equation 


dL_ 

Itt 


CiL 


L 

Cik2 


(A.5) 


Recalling that the volume of an EG is approximately 500 /xm^, closely packed 
cells have a density of 1 cell/500/xm^ = 2 x lO^cells/mm^. We assume the 
carrying capacity to be 10% of this, so that Cik 2 ~ 2 x lO^cells/mm^. In the 


following, we will fit experimental data to the solution of (A.5), which is 


L{t) = 


Cik2Loe^^^ 


cik2 + - 1 ) 


(A.6) 


^'Normal” proliferation rate Ci 

In [lOlj the different responses of lymphatic, venous and arterial endothe¬ 
lial cells to angiopoietins is studied; note that the cells were isolated and cul¬ 
tured from bovine mesenteric vessels. [1011 Figure 4B] shows the evolution 
of LEG density in time. After converting these data (in particular, those cor¬ 
responding to control, or 10% FBS) into cells/mm^, we can fit the function 
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(A. 6 ) to them, obtaining the estimate Ci = 0.42 day“^ with 95% confidence 
interval (0.15,0.70). Again, we can compare this value with other estimates 
obtained from different biological sources, such as [HaunniiHij; also, a similar 
parameter is estimated in [58]. Although the numerical values for ci found 
in these other references are all different, it is noticeable that they are all 
around 10“^. This make us very confident in estimating this parameter. 


Maximum density of cells (per unit time) k 2 

In the previous section we found different possible values for ci. Since 
Cik2 = 2 X lO^cells/mm^, we can easily obtain k2'. 

• From m we get k 2 ~ 4.71 x 10^ cells • day/mm^ (this is the value 
used in our model); 

• From [HH] we get /c 2 ~ 4.77 x 10^ cells • day/mm^ for bovine cornea 
and k 2 ~ 5.43 x 10^ cells • day/mm^ for bovine fetal heart; 

• From m we get k 2 ~ 3.13 x 10^ cells • day/mm^; 

• From [H7| we get k 2 ~ 4.48 x 10^ cells • day/mm^ for 10% FBS 
and k 2 ~ 18.52 x 10^ cells • day/mm^ for 2% FBS. 


VEGF-dependence of LECs growth C 2 ,C 3 

To estimate C 2 and C 3 , we consider only the exponential VEGF-dependent 
part of the LEG equation, that is 


dL _ f V \ 

dt \ ^ ^ C 2 + CsV) 


whose solution is 


L[t) — Lo • exp 



C2 + C3E J 



(A.7) 


(A. 8 ) 


Recall that we already have an estimate for Ci. Also, notice that, if time is 
fixed, (A. 8 ) can be seen as a function of V only. 

m Eigure 7A] shows the response of rat mesenteric LEGs to VEGE-A 
and VEGE-G at low serum conditions ( 2 % EBS). The cells were seeded at 
the density of 16,000 per well in 24-well plates. Then, VEGEies and mature 
VEGE-G (2-100 ng/mL) were added 4 hours after seeding. Einally, cells were 
counted 72 hours later. This provides a set of data giving the cell densities for 
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different concentrations of VEGF-A and C. These data refer to time t = 72 
hours = 3 days, and our approach is to fit the function (A. 8 ) as a function of 
V (with t fixed at 3 days) to the experimental data. Recalling that 1 ng/mL 
= 1 pg/mm^, and considering a standard well of 1 mL = lO^mm^ (for a 
24-well cell culture plate), we can convert these data into suitable units and 
use the MatLab function nlinfit to fit (A. 8 ) to them. This gives ^ = 42 
days and 63 = 4.1 pg/day/mm^, with 95% confidence intervals (—7.7, 92) and 
(2.9, 5.3) respectively. 


TGF(3-dependence of LECs growth C 4 

We estimate C 4 using experimental data obtained in the absence of VEGF. 
Therefore we consider only the part of the LEG equation concerning TGF-/3 
regulation of cell growth: 


— = ( ) L 

dt y 1 T C 4 T J 


whose solution is 


(A.9) 


L{t) — Lo • exp 


Cl 


1 + C4T 



(A.IO) 


In [ 88 ] the inhibitory action of TGF-/3 on bovine endothelial cells is studied. 
[881 Figures 1(a) and 1(b)] demonstrate the growth in time of bovine cornea 
and fetal heart endothelial cells respectively, while [881 Figure 1(c)] shows 
the inhibition on cell growth by TGF-/3. We can thus use these figures in 
different ways: 


In 


Figures 1(a) and 1(b)] the amount of TGF-/3 is fixed, so we 


can consider (A.IO) as a function of t only. Recalling that previously 
we found the value Ci = 0.42 day“^ for [HHl Figure 1 (a )], an d Ci = 
0.37 day“^ for [HHl Figure 1 (b)], we can fit the function (A.IO) to the 


data corresponding to T = 2 ng/mL = 2 pg/mm^ to find the parameter 
C 4 . The MatLab functions nlinfit and nlparci give us the estimate 
Q = 0.16 and its 95% confidence interval (—0.25, 0.58) for cornea EGs, 
and Q = 0.22 with 95% confidence (—0.50,0.95) for fetal heart EGs. 

Another strategy is to use the data contained in [HHl Figure 1 (c)], fixing 


the time {t — A days) and considering (A.IO) as a function of T only. For 
bovine cornea EGs, in this case we have 64 = 0.32 with 95% confidence 
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interval (0.12,0.51). Similarly the data for fetal heart ECs gives 64 = 
1.3 with 95% confidence interval (—0.66,3.2). 


Looking at the confidence interval for each of the estimates found above, we 
argue that the most “reliable” values for C 4 are the first three: 64 = 0.16, 
Q = 0.22 and q = 0.32. Hence, in our model we chose to take the average 
of these numbers, that is 64 = 0.24 mm^/pg. 

To compare this number with a similar estimate found in another source, 
we consider [891 Figure 1], which also shows how cell growth is influenced 
by TGF-^. Here bovine retinal and aortic endothelial cells were plated at 
25 cells/cm^ and TGF^-1 was added at different proliferation stages. Gell 
numbers were determined 5 days after the addition of different concentrations 
of TGF/3-1. Gonsidering the data in iHa Figure IB] and assuming a dish of 
10 mm height, we can fit the function (A. 10) to the data with the Mat Lab 
function nlinfit, taking Lq ~ 10 cells/mm^ (estimated from [891 Figure 
lA]) and consequently Ci = (1/^) ln(L(t)/Lo) = 0.37 day“^, since t = 5 days. 
This gives the estimate 64 = 6 with 95% confidence interval (—2.5,14). 


Threshold levels L* 

In [ 31 ] it is observed that in a wound space LEGs begin to organize in a 
network-fashion after about 25 days. Hence, we take L* to be the number of 
LEGs present at day 25 during normal repair. In [311 Figure 2 E] we see that 
at that time there are about 80 cells. This value corresponds to a section 
of width 12 /im. In addition, through IMl Figure 2D], we can calculate the 
observed wound area, which results to be 5.6 x 10^ /juT. In this way we get a 
volume of 0.0067 mm^ which contains 80 cells, which corresponds to 1.2 x 10^ 
cells/mm^. We take L* to be 10^ cells/mm^. 

For C*, we assume that LEGs stop coming into the wound when capillaries 
reach a level which is not far from the final one, that we have estimated to 
be = 1.2 X 10^ cells/mm^. Therefore we take C* = 10^ cells/mm^. Note 
that our estimate for L* and are the same because in [311 Figure 2E] it 
happens that the cell numbers counted at day 25 and at day 60 are about 
the same. 


Appendix B. Calculation of the steady states 

For T we have immediately 




di 


(B.l) 
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since the exponential tends to zero as t ^ oc. It follows that there is one 
T-steady state for every M-steady state. 

For M the situation is more complicated. Writing down the equation ([^ 
at the equilibrium and rearranging the terms one finds 


— /3—+ (/3r2 — d2 — pC)M + sm + olt — 

ki 02 + 1^ 




Plugging in the expression for T found in (B.l), the equation (B.2) becomes 

_ [ 0 , 10 ^] 

and M G [0,10^] 


a polynomial in M of degree 10. Plotting this polynomial for M G 

The graphs are 


gives the graphs shown in Figure B.20 


strongly indicative that there is only one stable steady state for (at least) 
M < 10^ and this is around 2000. This finding is in agreement with the 
estimate — 1875 cells/mm^ from the biological literature. 

The equilibrium equation for V is 


Sv + rsM - daV - -fVL = 0 ^ , (b.3) 

So there is one ^-equilibrium for every M, L equilibrium. 

At this point, it is more convenient to look at the capillary equation 
first, and afterwards consider the LEG one. There are two different cases to 
consider in the C-equation: 

• if L + C < L*, then cr = 0 and dC/dt — 0 always; 

• if L + C > L*, then cr = 1 and the equilibrium equation becomes 

(5i + 52V)L = 0 , 

whose only solution is L = 0 (since V cannot be negative). 

Coming to the L-equation, it is necessary to consider a number of different 
cases, since there are two piecewise-defined functions involved: 


CASE 

CONDITIONS 

FNs VALUE 

I 

L + C > L* and C < C* 

(t{L, C) = I and f{C) = 1 - C/C* 

II 

L + C > L* and C > C* 

a{L, C) = I and /(C) = 0 

III 

L + C < L* and C < C* 

(j(L, C) = 0 and /(C) = 1 - C/C* 

IV 

L + C < L* and C > C* 

a{L, C) = 0 and /(C) = 0 
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Figure B.20: Plots of the M-expresson at the equilibrium. The different lines 
refer to different representative values of C. 


Note that case IV seems to be not realistic^ since the estimates for the thresh¬ 
olds are L* ^ C*. 


CASE I Since we must have L = 0 in order to have a steady state in the C- 
equation, the equilibrium L-equation reduces to 


which implies 



bsV^ \ 
64 + V V 



= 0 , 


Ceq ^ 
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CASE II Again, we must have L = 0 for the equilibrium in the fifth equation. 

This time the steady state L-equation is automatically satisfied and 
therefore any value of C corresponds to a steady state. 


CASE III In this case, rearranging the fourth equation for L one gets 

1/ \ / 1 \ M + C" 


-^L^ + 
^2 


Cl + 



(B.4) 


where V depends on L and this dependence is given by the expression 


(B.3). Notice that if 7 = 0 then V does not depend on L. Hence, if we 


set 7 = 0 the expression (B.4) becomes a simple quadratic equation for 


L. In order to study how the system changes as V depends on L, we 


gradually increase 7 and see how the roots change. The function (B.4) 
is therefore plotted numerically versus L and the resulting graphs are 


reported in Figure B.21, These show that increasing 7 does not have a 


significant effect on the roots of (B.4): the function is always concave 


for 0 < L < 4 X 10^ and intercepts the horizontal axis once. Thus for 
C < C* (which is the case we are studying) there is only one intersection 
for L > 0 at about L — 2 x 10^. 
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Y= 0.0001 



Figure B.21: Plots of the expression (B.4) as a function of L for different 
values of 7 (more precisely, 7 = 0,10“^, 10“^, 1,10^, 10^). In each picture the 
function is plotted for three different representative values of C, in particular 
0 , CV 2 and C*. 
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